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Abstract 



CZ2 

In this study, we investigate the structure of the polar alkali-Strontium diatomic molecules as 



possible candidates for the realization of samples of new species of ultracold polar molecules. Using 
a quantum chemistry approach based on Effective Core Potentials and Core Polarization Potentials, 



03 

5/5 ■ we model these systems as effective three valence electron systems, allowing for calculation of 

O , 

electronic properties with Full Configuration Interaction. The potential curve and the permanent 

dipole moment of the S + ground state are determined as functions of the internuclear distances 



for LiSr, NaSr, KSr, RbSr, and CsSr molecules. These molecules are found to exhibit a significant 



permanent dipole moment, though smaller than those of the alkali-Rb molecules. 
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I. INTRODUCTION 



The prospects of the realization of quantum degenerate gases composed of polar systems, 
i.e. exhibiting a permanent dipole moment, either magnetic or electric, are among the most 
fascinating ones in the context of researches on ultracold matter |l( > related to many-body 
physics, quantum information, or quantum simulators of solid state physics. Indeed, neu- 
tral particles usually weakly interact via van der Waals short-range potentials varying with 
their mutual distance R as R~ 6 , which has been a key condition for the observation of Bose- 



a. 



In contrast, 



Einstein condensation in weakly interacting gases of ultracold alkali atoms 
particles with a permanent dipole moment interact through a long-range dipole-dipole po- 
tential varying as i? -3 which dominates the van der Waals interaction. Ultracold gases of 
such particles are predicted to evolve in a entirely new regime of strong interactions with 
pronounced anisotropy. For instance, the Chromium atom has a large magnetic dipole mo- 
ment of about 6/is, and once condensed [5|, |6(, the manifestation of anisotropic interactions 
has been observed experimentally [7], |8( . Ground state polar molecules, having permanent 
dipole moment (in the body-fixed frame) ranging between a few tenths of a Debye like LiNa 
or KRb, up to several Debye like LiCs [9j are actually strongly interacting with larger forces 
than in the Chromium case, and are expected to be better candidates for studying strongly 
interacting degenerate gases. Up to now, several attempts to observe ultracold ground state 



po 



ar molecules have been successful with bialkali molecules, namely KRb 10, lli|], LiCs 



121, NaCs 



131 ]. RbCs |14|. However the achievement of quantum degeneracy of ultracold 



molecules is conditioned by the phase space density of molecules in a unique quantum state, 
preferably the absolute ground state, defined not only by their vibrational and rotational 
quantum numbers, but also by their hyperfine structure 15Ml71]. 

As the number of atomic species other than alkalis, which are laser-cooled and subse- 
quently brought to quantum degeneracy is continuously increasing, the formation of novel 
species of polar molecules will soon be investigated in the ultracold regime 18J. For in- 



stance, photoassociation (PA) of Yb2 19 



and YbRb 



reported, while quantum degeneracy in Ytterbium [22|, 



2l| molecules 



ras already been 



23J, Calcium 24|, and Strontium 



25l . |26| atomic gases have recently been obtained. In the present paper, we investigate 



theoretically the structure of diatomic molecules composed of one Strontium atom and one 
alkali atom (Li, Na, K, Rb, Cs), for the first time for most of these species. In the follow- 
ing, we will refer to them as ASr molecules. In particular, we determined the permanent 
dipole moment of their ground state in order to provide guidance to the experimentalists 
to choose the appropriate combination of atomic species. In contrast with alkali dimers, 
such species also have a magnetic moment in the ground state due to their three valence 
electrons, which therefore makes them suitable for their manipulation with both electric and 
magnetic external fields. In Section HI1 we first recall our methodology based on quantum 
chemistry computations involving Effective Core Potentials (ECP), reducing the problem 
to an effective three valence electron system, for which Full Configuration Interaction can 
be achieved. Molecular potential curves and permanent dipole moments for the LiSr, NaSr, 
KSr, RbSr, and CsSr molecules in their electronic ground state are presented in Sections 
IIHI and IIVI emphasizing on the dipolar character of these molecules compared to bialkali 
molecules. 

II. METHOD OF CALCULATION 

We use the Configuration Interaction by Perturbation of a mult i- configuration wave func- 
tion Selected Iteratively (CIPSI) set of programs developed at the Universite Paul Sabatier 
in Toulouse (France). The CIPSI package incorporates both the use of Effective Core Poten- 



tials (ECP)|27|,|28|, and ^-dependent Core Polarization Potentials (CPP) |29|,|30(. Effective 
core potentials model the effect of inner shells electrons into a semi-empirical potential thus 
reducing the number of active electrons on each atom, namely one for the alkali atom, and 
two for the Strontium atom. Core polarization potentials embody the polarization effects of 
these effective cores with the valence ones. Thanks to the reduced number of active electrons, 
it becomes possible to perform a Full Configuration Interaction (FCI) in the valence space. 
This gives us a number of excited molecular states in addition to the lowest state in each 
symmetry. This method has proven to be efficient and accurate for alkali diatomics which 
are treated as two-electrons system [9|, |3jJ, and we modified the CIPSI package accordingly 
for effective three-electron systems. 



Effective core potentials are from ref. 32] . For the strontium atom, we use the so- 
called "large core" ECP for the 36 electrons of the Sr 2+ core, and basis sets composed of 



uncontracted Gaussian functions. The basis sets for alkali atoms are taken from ref.[9(, where 
/ gaussian have been removed to reduce the size of the configuration space. We checked on 
the alkali atoms that this has no influence on the lowest energy levels we are interested in 
here. The basis set for the strontium atom is taken from [33]. We actually used two slightly 
different basis sets: we first designed sets for ASr + ions, i. e. 8s7p4d for Lithium, 7s6p5d for 
Sodium, 7s5p7d for Potassium, 7s4p5d for Rubidium, 7s4p5d for Cesium, and 7s6p7d for 
Strontium ion whose exponents are reported in Table fl] These results will be discussed in a 
separate publication. For the neutral systems ASr, we reduced the Sr + set by a few orbitals 
down to 5s5pQd, due to size constraints. The influence on the atomic energy levels remains 
small as it can be seen below. 

In total, the number of configurations for the diatomic calculations is larger than 10 5 , and 
depends on the molecular symmetries. We use a recently modified version of the ^-dependent 
core polarization potential initially proposed by Foucrault et al. [30(] , which yields the correct 



asymptotic behavior of the potential energy at large internuclear distance 34J , in contrast 
with most of the previous calculations performed on alkali dimers with this method. A set 
of cutoff radii, reported in Table [Til are optimized to reproduce experimental values for the 
lowest s, p, and d atomic levels of the alkali atoms and of the strontium cation (Table ITTTI) . 



Polarizabilities for the relevant cores Li + , Na + , K + , R 



d + , Cs + ,and Sr 2+ cores are taken from 



experimental and semi-empirical values of refs. |35l. [36|. 

Table IIHI shows that such an optimization yields a satisfactory representation of the 
Strontium cation spectrum. The corresponding fitted levels of the alkali atoms are not 
reported here, as we are mainly interested into the lowest level of s, p, and d symmetry 
whose energies are represented exactly after the tuning of the cut-off radii of Table HTl 
Therefore the energies of a pair composed by one alkali atom and a Strontium cation in 
their lowest levels at infinite separations are well modeled here. 

III. THE GROUND STATE POTENTIAL CURVE OF ALKALI-STRONTIUM 
MOLECULES 

As dissociation limits of ASr molecules involve a neutral alkali atom and a neutral Stron- 
tium atom, we first performed a first FCI calculation with the small and large basis sets and 



Li 



Na 



K 



Rb 



Cs 



Sr^ 



p s 1.4519565 1.450275 2.10816 2.5365 2.85585 2.13999 
p p 1.013 1.6452591 1.9722 2.279 2.544 2.183 
Pd 0.6 1.5 1.9405 2.5217 2.89815 1.706 



Lr 



Na+ 



K" 



Rb H 



Cs^ 



Sr 



2+ 



a 0.1997 0.9987 5.472 9.245 16.33 5.67 



TABLE II. Optimized cutoff radii p^ (in atomic units ao) and core polarizabilities a of the relevant 



ionic cores Li+, Na+, K+, Rb+, Cs+ [35], Sr 2+ 



361 ] entering into the Core Polarization Potentials 



(CPP). 



angular momentum 


Li 


Na 


K 


Rb 


Cs 


Sr+ 


s 


2.464 


2.8357 


0.9312 


1.292561 


0.328926 


0.79174 




1.991 


0.49318 


0.2676 


0.824992 


0.241529 


0.316178 




0.582 


0.072085 


0.0427 


0.234682 


0.050502 


0.066565 




0.200 


0.036061 


0.03915 


0.032072 


0.029302 


0.02699 




0.070 


0.016674 


0.01448 


0.013962 


0.013282 


0.013495 




0.031 


0.00693 


0.0055 


0.00575 


0.00528 


0.009 




0.015 


0.00287 


0.0026 


0.0025 


0.003 


0.003 




0.007 












P 


0.630 


0.431 


0.133 


0.128 


0.1 


0.225825 




0.240 


0.09276 


0.05128 


0.040097 


0.0405 


0.095691 




0.098 


0.03562 


0.01642 


0.014261 


0.0162 


0.042077 




0.043 


0.01447 


0.0052 


0.00485 


0.00443 


0.018077 




0.020 


0.0058 


0.0022 






0.0090385 




0.010 


0.00023 








0.005 




0.005 












d 


0.180 


0.292 


1.255 


0.408807 


0.196894 


3.618081 




0.080 


0.06361 


0.343 


0.096036 


0.067471 


0.999656 




0.22796 


0.02273 


0.109 


0.026807 


0.027948 


0.390735 




0.008574 


0.008852 


0.0294 


0.009551 


0.010712 


0.12277 






0.00352 


0.01013 
0.0039 
0.0018 


0.004 


0.003 


0.036655 

0.018327 

0.009 


symmetry 


LiSr 


NaSr 


KSr 


RbSr 


CsSr 




l s+ i A x 2_ y 2} 


129945 


132792 


161703 


110337 


110337 




n 


120032 


123684 


151478 


102540 


102540 




{^ ,AxyS 


110240 


114840 


142138 


95457 


95457 





TABLE I. Gaussian basis set exponents used to represent the effective one-electron systems of 
the present work, and resulting number of configurations representing the effective three-electron 
diatomic systems A-Sr (with A = Li, Na, K, Rb, Cs) for the relevant molecular symmetries. The 
actual basis sets used for ASr neutral molecules include all exponents but the ones in italic, while 
the whole set has been used for ASr + systems. The three lowest lines provide the size of the 
hamiltonian for the Full Configuration Interaction. 



State 


Exp [39] 


(a) 


A E 


(b) 


A £ 


[40] 


A E 


[41] 


A E 


5s 


-0.4053569 


-0.405354 





-0.405358 





-0.405522 


36 


-0.40839 


666 


4d 


-0.3382685 


-0.338269 





-0.338269 





-0.338247 


-4 


-0.341908 


799 


5p 


-0.2948674 


-0.294868 





-0.294869 





-0.294854 


-2 


-0.29588 


222 


6s 


-0.1878519 


-0.187521 


-73 


-0.187544 


-68 


-0.187469 


-84 






5d 


-0.1623282 


-0.160856 


-323 


-0.160947 


-303 


-0.161237 


-239 






6p 


-0.1503744 


-0.150234 


-31 


-0.150237 


-30 


-0.150150 


-49 






7s 


-0.1093568 






-0.109138 


-48 










6d 


-0.9758806 






-0.096543 


-230 










7 P 


-0.921336 






-0.091920 


-46 










8s 


-0.7166203 






-0.075047 


-34 











TABLE III. Binding energies (in a.u.) of Sr + levels obtained with (a) the small basis set, and 



(b) the large basis set for Sr + (Table [TJ> , and their deviations Ae (in cm x ) from experiment 



The ionization limit of Sr + is taken at 88965.18 cm 1 



39]. 



421 ] to calculate the experimental binding 



energies in a.u.. We use the mass-corrected Rydberg constant 109736.63 cm . The theoretical 
values of 4l( have been spin-averaged. 



involving two electrons, in order to calculate the lowest strontium energy levels. Results are 
given in Table IIV} and are compared to experimental values. The differences range between 
less than 200 cm -1 for the ground state up to about 450 cm -1 for the lowest excited levels, 

n 

as it was the case also for the Magnesium atom [37J. We note that such an accuracy is 
typical for this kind of calculations, as illustrated also by ECP+CPP FCI calculations on 
Calcium [38| . The use of the small basis set does not significantly affect the quality of the 

results. 

We display in Figure [T|a) the potential energy curves (PEC) for the AT 2 £ + ground state 
of the ASr systems, relative to the same origin of energies fixed at the dissociation limit 
Ans+Sr(5s 2 1 S) for comparison purpose. The corresponding main spectroscopic constants 
are reported in Table |V] The Strontium atom in its ground state is a closed-shell atom, 
so that the alkali-strontium molecules are expected to have a van der Waals character with 
a somewhat large equilibrium distance and weak binding energy, compared for instance 
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Exp 


(a) 


(b) 


Exp 


(a) 


(b) 




rpexp 


E b 


Sb 


E b 


4 


T^exp 


E e 


Se 


E e 


Se 


5s 2 x 5 


-0.61464148 


-0.615471 


-182 


-0.615484 


-184 

















5s5p 3 P° 


-0.54765155 


-0.548880 


-269 


-0.548899 


-273 


14702.6 


14613 


-87 


14613 


-88 


5sM 3 D e 


-0.53147099 


-0.529868 


351 


-0.529907 


343 


18253.7 


18779 


533 


18781 


528 


5s4d 1 D e 


-0.52283208 


-0.520684 


471 


-0.520733 


460 


20149.7 


20793 


653 


20795 


645 


5s5p 1 P° 


-0.51577530 


-0.517756 


-434 


-0.517792 


-442 


21698.5 


21447 


-252 


21440 


-257 



TABLE IV. Computed Strontium binding energies E b (in a.u. relative to the double ionization 
threshold of Strontium) and excitation energies E e (in cm -1 ) compared to the experimental values 
of ref. [39] (E^ xp and E^ xp ), and their differences (<5& = E b — E^ xp and 5 e = E e — E% xp ) in cm -1 , 
resulting from a two-electron FCI calculation with (a) the large basis set, and (b) the small basis 
set for Sr + (Table UJ). The experimental values are averaged over fine structure energies when 
appropriate. We use the mass-corrected Rydberg constant 109736.63 cm . 



to alkali diatomics with comparable masses like alkali-rubidium molecules. This feature is 
clearly visible in the figure, as the typical binding energies are three to five times smaller 
than for the related A-Rb molecules. The equilibrium distances are also systematically larger 
for ASr than for ARb molecules, except for LiSr. The LiSr molecule is a bit peculiar, as its 
ground state is quite deep, with an equilibrium distance comparable to the LiRb one. 



The potential depths of the ASr molecules immediately show that in ultracold conditions 

;hese molecules will not be stable against collisions between them, just like KRb molecules 

43[. The reaction channel yielding a Sr2 molecule and an alkali dimer is energetically open 

for all alkali species. Indeed, the Sr 2 well depth is about 1061 cm -1 44[, while those of the 



alkali dimers are 8517 cm 1 |45l, 6022 cm 1 



46|, 4451 cm" 1 J47J, 3993 cm" 1 J48j, 3649 cm" 1 



49j, for Li 2 , Na 2 , K 2 , Rb 2 , Cs 2 , respectively. 



From the long range part of these curves we can estimate the coefficient Cq of the leading 
term of the van der Waals interaction 




Internuclear distance (in units of a ) 



FIG. 1. (a) Potential energy curves and (b) permanent dipole moments (PDM) for the A" 2 S + 
ground state of the ASr molecules (A=Li, Na, K, Rb, Sr). Vertical lines guide the eye to locate 
the value of the PDM at the equilibrium distance. 

IV. THE PERMANENT DIPOLE MOMENT OF ALKALI-STRONTIUM MOLECULES 



As stated in the introduction, an important property of these molecules, in the perspective 
of cooling and trapping them at ultracold temperatures is the permanent dipole moment 
(PDM), which are drawn in Figure Hfb) and reported as well in Table |V] The PDM for 
ground state ASr molecules is predicted to be noticeable, even if they do not reach such 
high values than those of the alkali dimers [9|. It is worthwhile to mention that the relation 
between the magnitude of the PDM and the difference in the mass (and then in the size) of 
the constituting atoms is inverted compared to alkali dimers: the PDM of LiCs and LiRb 
are the largest ones among all alkali pairs, while the PDM of LiSr is by far the smallest 
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D e (cm- 1 ) 


R e (a.u.) 


LO e (cm -1 ) 


B e (cm- 1 ) 


d e (Debye) 


d (Debye) 


LiSr 


2587 


6.57 


184.9 


0.21 


0.34 


0.34 


NaSr 


1597 


7.22 


85.4 


0.063 


0.62 


0.63 


KSr 


1166 


8.33 


52.4 


0.032 


1.52 


1.50 


RbSr 


1073 


8.69 


32.3 


0.018 


1.54 


1.53 


CsSr 


1084 


9.06 


33.8 


0.013 


1.91 


1.91 


LiRb 


13759 


6.49 


210 




4.42 




NaRb 


5076 


6.82 


107 




3.41 




KRb 


4199 


7.63 


76 




0.62 




CsRb 


3907 


8.29 




1.28 







TABLE V. Main spectroscopic constants for the lowest electronic states of the ASr molecules 
(A=Li, Na, K, Rb, Sr): well depth D e , equilibrium distance R e , harmonic constant oj e , and 
rotational constant B e . The value of the dipole moment d e at the equilibrium distance, and averaged 
for the v = level (do) are also reported. The well depth, equilibrium distance, and permanent 
dipole moment at equilibrium of the related ARb molecules are displayed for comparison, following 



our computations reported in Refs. 



y,y, 



50] 



one compared to the other ASr species. This feature is amplified by the mismatch of the 
minimum of the PEC of LiSr and NaSr, and the maximum of the PDM curve, in contrast 
with the other species. 

V. CONCLUSION 



In this paper we investigate for the first time the electronic structure of the alkali- 
Strontium molecules with a quantum chemistry approach which accurately treat correlation 
between valence electron through a Full Configuration Interaction, while the correlation be- 
tween valence and core electrons is accounted for through effective potentials. This approach 
has proven to be powerful for similar calculations on alkali dimers. The potential curve and 
the main spectroscopic constants, as well as the permanent dipole moment (PDM) of the 
2 S + ground state have been determined for LiSr, NaSr, KSr, RbSr, and CsSr molecules. 
Though their PDM is generally smaller than the alkali-Rb molecules, these systems are 

11 



possible candidates for achieving a new kind of ultracold molecular sample with anisotropic 
interactions. Just like heteronuclear alkali dimers, ultracold molecules in their ground state 
can probably be created by photoassociation, and this will be the topic of a further study. 
While completing this paper, Zuchowski, Aldegunde and Hutson informed us about their 



study about the structure of the RbSr ground state [5l|], using an all-electron quantum 
chemistry approach. Their results are in good agreement with ours. For instance, they ob- 
tained an equilibrium distance R e = 8.86 a.u., a well depth of about 1000 cm" 1 , and a PDM 
of 1.36 D at R e (compared to our values R e = 8.69 a.u., D e = 1073 cm -1 , and d e = 1.54 D). 
The slight differences are probably due to the sensitivity of the results on the treatment of 
electron correlation, and further experiment should help to clarify the situation. In their 
paper, Zuchowski et al. demonstrated the existence of Feshbach resonances in RbSr. There- 
fore, our next study on excited states of these molecules should guide experimentalists to 
design ways for achieving adiabatic population transfer from molecules created by Feshbach 
association down to the lowest energy level of the molecule, just like it has been recently 



done for Cs2 



17 



52|and KRb 



15[ molecules. 
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